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ABSTRACT 

Two-point  correlations,  energy  spectra,  and 
length  scales  are  examined  in  the  vicinity  of  a  free 
surface,  modeled  as  a  shear-free  boundary  in  a  di¬ 
rect  numerical  simulation  of  open  channel  flow.  The 
length  scale  results  indicate  that  a  typical  eddy  is  flat¬ 
tened  as  it  interacts  with  the  surface.  The  scales  as¬ 
sociated  with  the  vertical  component  of  velocity  seem 
to  determine  the  extent  of  the  source  layer  described 
in  the  Hunt-Graham  model.  The  energy  spectra  show 
qualitative  agreement  with  the  model,  though  higher 
resolution  calculations  will  be  required  to  make  more 
quantitative  comparisons.  Additionally,  the  proxim¬ 
ity  of  the  free  surface  to  the  bottom  solid  wall  of  the 
channel  evidences  itself  as  a  wall-layer  streaky  struc¬ 
ture  which  persists  to  a  noticeably  greater  distance 
away  from  the  wall.  Some  speculations  are  offered  to 
explain  this  effect.  ^ 

NOMENCLATURE 

=  channel  height  ^ 

=  turbulent  kinetic  energy 
=  u/ut,  viscous  length  scale 
=  U,h/v,  Reynolds  number 
=  xifhlv,  wall  Reynolds  number 
=  Uoo^lv,  momentum  thickness 
Reynolds  number 
=  I'/uJ,  viscous  timescale 
=  instantaneous  velocity  vector 
=  instantaneous  velocity  component 
=  fluctuating  velocity  component 
=  friction  velocity 

=  turbulent  microscale 
=  turbulent  macroscale 
=  two-point  correlation  tensor 
=  coordinate  direction 
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(  ) 

—  averaged  quantity 

=  specral  density 

=  streamwise  wavenumber 

U 

=  kinematic  viscosity 

P 

=  density 

e 

=  momentum  thickness 

Tw 

=  shear  stress 

n 

iii 

=  instantaneous  vorticity  vector 
=  instantaneous  vorticity  component 

Subscripts 

i 

s 

oo 

w 


—  1,2,3,  coordinate  directions 
=  value  at  free  surface 
=  value  in  free  stream 
=  value  at  wall 


1.  INTRODUCTION 


The  study  of  the  structure  of  turbulence  near  a 
free  surface  is  obviously  important  to  our  understand¬ 
ing  of  the  complex  interaction  of  the  atmosphere  and 
upper  ocean.  It  is  also  of  fundamental  relevance  to  the 
wall-bounded  turbulence  problem,  since  it  isolates  the 
boundary  influence  on  turbulent  fluctuations  from  the 
turbulence  production  mechanism  at  the  wall.  The 
first  detailed  experiment  which  addressed  itself  to  this 
particular  problem  was  that  of  Uzkan  and  Reynolds* 
(UR).  They  passed  grid  generated  homogeneous  tur¬ 
bulence  over  a  wall  which  moved  with  the  mean  flow 
and  therefore  generated  no  mean  shear  at  the  bound¬ 
ary.  They  found  that  the  streamwise  turbulence  in¬ 
tensity  near  the  shear-free  boundary  did  not  peak  as 
it  does  near  a  stationary  solid  wall,  but  instead  de¬ 
creased  monotonically  from  its  free  stream  value  to 
zero  at  the  boundary.  Later,  Thomas  and  Hancock^ 
(TH)  performed  a  similar  experiment  at  a  Reynolds 
number  about  20  times  greater  than  that  of  UR  and 
found  that  the  intensity  of  the  streamwise  component 
increases  as  the  boundary  is  approached. 

The  discrepancy  between  these  two  results  was 
explained  satisfactorily  by  Hunt  and  Graham^  (HG) 
who  proposed  a  two  layer  model  for  the  interaction. 
At  high  turbulent  Reynolds  numbers  there  exists  a 
thin  viscous  layer  near  the  wall  embedded  in  a  larger 
source  layer.  The  source  layer  should  be  roughly  the 
size  of  the  integral  length  scale  of  the  free  stream  tur¬ 
bulence  and  exists  essentially  because  of  the  no  mass- 
flux  condition  at  the  boundary.  Their  theory  predicts 
a  redistribution  of  turbulent  energy  in  the  source  layer 


from  the  vertical  component  of  velocity  to  the  stream- 
wise  and  spanwise  components.  The  UR  result  was 
easily  explained  since,  at  the  low  Reynolds  number  of 
their  experiment,  the  viscous  layer  dominated  the  re¬ 
gion  near  the  surface  and  the  turbulence  was  accord¬ 
ingly  damped.  At  the  higher  Reynolds  number  of  the 
TH  experiment,  the  source  layer  dominated  and  the 
redistribution  of  the  turbulence  behaved  according  to 
the  HG  model.  Many  of  these  results  were  later  con¬ 
firmed  by  the  large  eddy  simulations  of  Biringen  and 
Reynolds'*.  Recently,  Brumley  and  Jirka®  (BJ)  pre¬ 
sented  results  for  experiments  in  which  homogeneous 
turbulence  interacted  with  a  free  surface.  Their  re¬ 
sults  agreed  reasonably  well  with  a  modified  form  of 
the  HG  model. 

The  simulations  performed  here  were  designed  to 
represent  as  closely  as  possible  the  physics  of  free  sur¬ 
face/turbulence  interaction  in  which  the  effects  of  sur¬ 
face  waves  can  be  safely  neglected.  For  this  purpose, 
fully  developed  turbulence  between  a  solid  wail  and  a 
free  surface  is  simulated.  The  physical  processes  rep¬ 
resented  by  these  simulations  differ  in  some  important 
respects  from  processes  involved  in  the  physical  exper¬ 
iments  noted  above.  First,  in  these  simulations,  no 
viscous  layer  can  develop  since  ut  and  «3,  the  fluctu¬ 
ating  streamwise  and  spanwise  velocity  components, 
are  not  forced  to  zero  as  in  the  UR  and  TH  experi¬ 
ments.  In  fact,  even  in  the  BJ  experiments  a  viscous 
layer  developed  near  the  surface  due  to  the  presence 
of  surface  contaminants.  Secondly,  in  these  simula¬ 
tions,  the  turbulence  impinging  on  the  free  surface  is 
not  isotropic  since  it  is  being  generated  at  a  solid  wall. 
It  is  evident  that  the  solid  wall  is  acting  as  a  genera¬ 
tor  of  anisotropic  turbulence  which  is  then  converted 
toward  the  free  surface.  Lastly,  it  is  possible  to  ex¬ 
amine  the  question  of  the  influence,  if  any,  of  the  free 
surface  on  the  turbulent  structure  at  the  solid  surface 
boundary.  In  this  work  two-point  correlations,  energy 
spectra,  and  turbulent  length  scales  will  be  examined 
in  an  effort  to  investigate  the  turbulent  structure  near 
a  free  surface  in  the  absence  of  surface  waves. 

2.  DIRECT  NUMERICAL  SIMULATION 


The  incompressible  three-dimensional  Navier- 
Stokes  equations  are  solved  for  initial  and  boundary 
conditions  approximating  a  turbulent  open  channel 
flow  of  water  at  a  Reynolds  number,  Re^,  based  on 
the  channel  height,  h,  and  the  mean  velocity  at  the 
free  surface,  U,,  of  2340.  The  notation  Zi,  X2.  and 
£3  is  used  to  denote  the  streamwise,  wall-normal, 
and  spanwise  coordinates  respectively.  The  govern¬ 
ing  equations,  formulated  in  the  manner  suggested 
by  Orszag  and  Patera^  and  later  implemented  in  a 
.simpler  form  by  Kim,  Mom  and  Moser’,  con.sist  of  a 


fourth  order  equation  for  the  vertical  velocity,  C/j  . 
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and  a  second  order  equation  for  the  vertical  vortic- 
ity,  122.' 


(s  -  S  «>  =  a?;'^  X  ^ 

(2) 

where  all  variables  are  made  non-dimensional  by  h 
and  the  initial  value  of  U,  .  Here,  the  instantaneous 
velocity  vector  is  given  by  U  and  the  instantaneous 
vorticity  vector  is  defined  by  fl  =  (V  x  U).  Follow¬ 
ing  the  solution  of  equations  1  and  2,  the  streamwise 
and  spanwise  velocity  components,  Uj  and  U3.  are 
recovered  from  the  incompressibility  condition. 

The  equations  of  motion  are  solved  in  Fourier- 
Chebyshev  space  where  Fourier  modes  are  employed 
in  the  horizontal  plane  and  Chebyshev  modes  in  the 
wall  normal  direction.  The  calculations  are  performed 
on  a  64  X  65  X  48  grid  in  xi,Z2,Z3  respectively. 
With  the  geometry  scaled  by  the  channel  height,  the 
streamwise,  vertical,  and  transverse  dimensions  of  the 
channel  are  4ir,  I,  and  3ir/2  respectively.  In  terms  of 
the  viscous  parameters  consisting  of  the  friction  ve¬ 
locity,  Ur.  and  the  kinematic  viscosity,  i/,  the  domain 
is  1684£*  X  134f*  x  632f*  and  the  Reynolds  number, 
/?*,  is  134.  To  facilitate  substantive  qualitative  com¬ 
parisons  with  the  wall-bounded  turbulence  problem, 
a  companion  calculation  for  a  closed  channel  flow  is 
utilized.  For  reasons  of  economy,  this  calculation  is 
at  half  the  wall-normal  resolution  of  the  open  chan¬ 
nel  case  and  is  at  a  somewhat  lower  Reynolds  number, 
R'  =  125.  Nevertheless  the  behavior  of  these  data  for 
all  aspects  examined  is  identical  to  that  reported  by 
KMM. 

The  boundary  conditions  utilized  are  periodic  on 
all  dependent  variables  in  the  streamwise  and  span- 
wise  directions.  No-slip  conditions  are  used  at  the 
channel  bottom  while  the  free  surface  is  approximated 
as  a  rigid,  free-slip  surface  with  vanishing  shear.  The 
shear-free  rigid  lid  condition  is  an  approximation  to 
the  exact  free  surface  condition  which  is  valid  at  low 
Froude  number  for  a  surface  free  of  any  contaminants. 
Leighton  et  al.^  have  estimated  the  surface  displace¬ 
ment  a  posteriori  from  the  results  of  the  simulation 
using  the  channel  height,  h  =  0.04m,  from  the  experi¬ 
ments  of  Komori  et  ai.^  .  For  these  conditions  the  rm.s 
surface  deflections  are  expected  to  be  approximately 
16  X  lO'^m  (0. 004/1)  and  are  clearly  negligible  as  ob¬ 
served  in  the  laboratory  experiments  The  boundary 
conditions  at  the  solid  wall  {x^  =0)  and  the  free  sur¬ 
face  1x2  =  1 )  are. 
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U\  —  =  Ua  —  0;  f  2  =  0, 


and 


dLh 

dX2 


dUz 

dxo 


=  t/2  =  0;i2  =  1. 


(3) 

(4) 


The  code  developed  for  the  simulation  is  de¬ 
signed  to  run  on  the  CRAY  X-MP/24  computer  at 
the  Naval  Research  Laboratory  and  requires  approxi¬ 
mately  10~®  CPU  seconds  per  timestep  per  grid  point. 
After  the  wall  shear  stress  achieved  a  statistically 
steady  behavior,  42  realizations  of  the  instantaneous 
velocity  data  were  saved  during  a  time  interval  of  ap¬ 
proximately  4000t*  where  t‘  =  v/u?.  Statistics  were 
obtained  by  averaging  in  the  streamwise  and  spanwise 
directions  as  well  as  over  all  realizations.  Swean  et 
a/.^°  find  good  agreement  between  these  simulations 
and  open  channel  flow  experiments. 

3.  TWO-POINT  CORRELATIONS 
AND  ENERGY  SPECTRA 

The  turbulent  structure  near  the  free  surface  is 
revealed  in  some  detail  by  examining  the  two-point 
correlations  and  energy  spectra  at  different  depths  be¬ 
low  the  surface.  The  two-point  correlation  function, 
Rij,  is  defined  by  : 

Ajr3,X2,X2')  = 

U,(li,X3.  J2)U;(Xi',X3',X2') 

where,  x/  =  xy  -f-  Axy,  j=l,3.  Here,  only  the  prop¬ 
erties  of  Rij  for  which  X2  =  x^  and  i  =  j  will  be 
described.  These  correlations  were  computed  by  av¬ 
eraging  over  all  flow  realizations  and  all  flow  symme¬ 
tries  (see  Sirovich^’).  Figure  1  shows  a  comparison  of 
the  streamwise  correlations  for  open  and  closed  chan¬ 
nel  flows  at  several  X2  locations.  Note  that  in  each 
figure  a  secondary  axis  is  given  showing  the  correla¬ 
tion  length  in  terms  of  wall  variables,  i.e.  Ax^  for 
this  case.  In  the  region  close  to  the  bottom  solid  wall 
(0  <  X2  <  0.5;  0  <  Xj  <  68)  the  streamwise  correla¬ 
tions  are  virtually  identical  for  all  three  velocity  com¬ 
ponents,  As  an  example,  note  the  similarity  between 
the  correlations  at  X2  =  0.071  for  open  channel  flow 
with  those  of  closed  channel  flow  at  X2  =  0  076.  The 
only  notable  difference  is  a  somewhat  longer  stream- 
wise  correlation  length  for  Ui  in  the  open  channel  case. 
At  di.stances  farther  from  the  wall,  however,  the  dif¬ 
ferences  between  the  two  flows  become  increasingly 
more  pronounced.  As  the  free  surface  is  approached, 
two  trends  are  evident.  First,  the  streamwise  dis¬ 
tance  at  which  /?33( Axi ,  0,  X2),  (subsequently  de¬ 
noted  /?33(Axi))  attains  its  minimum  value  increases 
as  the  free  surface  is  approached.  For  example,  at 
X2  =  0  524.  the  minimum  occurs  at  Axi  ~  1.34  and 


as  the  free  surface  is  approached  this  increases  roughly 
17%  to  Axi  ~  1.57.  For  the  closed  channel  case,  how¬ 
ever,  this  distance  undergoes  virtually  no  change  as 
X2  varies  from  0  572  to  the  centerline.  Secondly,  we 
observe  that  the  streamwise  correlation  length  of  the 
vertical  component  of  velocity  in  the  open  channel 
case,  which  may  be  loosely  defined  by  the  first  zero 
crossing,  decreases  significantly  between  xj  =  0.524, 
where  no  zero  crossing  exists,  and  X2  =  0.952  where  it 
attains  a  value  of  Axi  ~  0  75  .  In  the  closed  channel 
case  this  scale  also  decreases,  but  not  nearly  as  rapidly 
as  in  open  channel  flow  where  there  appears  to  be  a 
clearer  separation  between  /?22(-^^i)  and  i?33(Axi) 
near  the  free  surface.  In  Fig.  2  the  corresponding 
results  for  the  spanwise  correlations,  /Zyj(Ax3),  are 
shown  for  both  cases.  As  in  the  streamwise  results, 
close  to  the  wall  there  appear  to  be  no  significant  dif¬ 
ferences  between  these  flows.  Farther  from  the  wall 
(x2  >  0.5),  it  is  evident  that  the  correlation  length 
for  the  U2  component  in  the  open  channel  flow  is  sig¬ 
nificantly  smaller  than  those  for  the  other  two  veloc¬ 
ity  components.  This  behavior  is  not  evident  for  the 
closed  channel  case  in  this  same  region.  Again,  as 
with  the  streamwise  correlation  length  results,  there 
is  a  clear  separation  between  /?22(  Axa)  and  Rzs^Axz) 
as  the  free  surface  is  approached.  One  feature  of  note 
is  the  existence  of  a  discernable  local  minimum  in  the 
spanwise  correlation  of  the  streamwise  velocity  com¬ 
ponent  out  to  normal  locations  as  large  as  12  0  8 

(i.e.,  ss  108  viscous  lengths)  in  the  open  channel  case. 
This  indicates  a  considerable  persistence  of  a  spanwise 
periodic  structure  in  the  flow  which  will  be  discussed 
in  more  detail  in  Section  5 

Energy  spectra,  ^ij,  as  a  function  of  the  stream- 
wise  wavenumber,  «),  are  shown  in  Fig.  3  for  two 
different  surface  normal  locations  in  close  proximity 
to  the  free  surface  boundary  (x2  =  0.978  and  0  952). 
Additionally,  a  third  energy  spectrum  at  xt  =  0.798 
is  given  in  each  plot  as  a  reference  condition.  Ex¬ 
amination  of  the  turbulence  intensity  profiles  (refer 
to  the  Fig.  7  discussed  later)  at  this  location  indi¬ 
cates  that  free  surface  effects  will  be  negligible.  The 
‘^n(ai)  spectra  show  that  very  near  the  free  surface, 
the  energy  at  low  wavenumbers  remains  unchanged. 
However,  for  the  intermediate  band  (2  <  xi  <  10)  a 
small  incre^ise  in  energy  is  evident.  There  is  no  change 
in  the  spectra  at  high  wavenumbers.  These  features 
are  quite  consistent  with  the  HG  predictions  From 
the  <^22(^1 )  spectra  it  is  quite  evident  that,  as  the  free 
surface  is  approached,  the  energy  at  low  wavenumbers 
decreases  more  rapidly  than  at  high  wavenumbers. 
This  is  also  consistent  with  HG.  However,  it  is  to  be 
noted  that  the  HG  model  predicts  that  the  '^22(''i) 
spectrum  far  from  the  boundary  will  merge  with  the 
spectra  near  the  free  surface  at  a  wavenumber  of  or¬ 
der  2x7 Ax2,  where  Ax2  is  the  vertical  distance  from 
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the  free  surface.  This  is  effectively  due  to  the  cut¬ 
ting  off  of  eddies  smaller  than  Axj  by  the  presence  of 
the  surface.  In  these  calculations,  however,  the  min¬ 
imum  resolvable  streamwise  length  scale  dictated  by 
the  grid  spacing  is  4jr/64  =  0.196,  which  is  in  fact  of 
the  same  order  as  the  source  layer.  It  is  evident  that 
higher  resolution  will  be  required  to  resolve  these  ef¬ 
fects.  Nevertheless,  the  basic  structure  of  is 

suggestive  of  this  kind  of  wavenumber  cutoff  behav  ior. 
The  ^33{ki)  spectra  show  some  increase  in  energy  at 
low  wavenumbers  as  the  free  surface  is  approached, 
but  virtually  no  change  for  fci  >  3.  These  observa¬ 
tions  are  also  qualitatively  consistent  with  HG. 

4.  MACRO  AND  MICROSCALES 

Turbulent  flows  are  known  to  contain  a  wide 
range  of  length  scales;  here  we  examine  both  the  tur¬ 
bulent  macroscales  and  microscales.  The  macroscale 
can  be  considered  the  length  scale  that  represents 
the  size  of  a  typical  energy  containing  eddy  which 
is  eventually  broken  up  and  dissipated  by  viscosity 
at  smaller  scales^^.  The  microscale,  though  not  the 
smallest  length  scale  in  the  flow,  can  be  thought  of 
as  an  average  length  within  which  most  of  the  energy 
dissipation  occurs.  In  high  Reynolds  number  flows 
there  is  a  large  separation  between  these  two  scales, 
but  in  the  current  computations  this  separation  is  not 
large. 

The  microscale,  A,  corresponding  to  velocity 
component  Uj  in  direction  x,  is  defined  by; 
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(6) 


If  the  turbulence  is  homogeneous  in  direction  Xj  then 
it  can  be  shown  that  an  equivalent  definition  is  : 


The  xi  and  X3  microscales  are  computed  using  both 
definitions  given  above  and  nominally  produce  iden¬ 
tical  results  The  macroscale  A  is  defined  by  : 


the  largest  eddy  structures  near  the  wall.  For  these 
reasons,  the  microscales  reported  below  may  actually 
be  larger  in  some  cases  than  the  macroscales.  This  is 
simply  an  artifact  of  the  method  used  here.  A  better 
estimate  for  the  macroscales  can  be  had  in  some  cases 
by  using  the  length  associated  with  the  second  zero 
crossing  of  the  correlation  function. 

Figure  4  shows  the  results  of  the  calculation  of 
the  streamwise  macroscale,  A 1  (i.eAi^),  and  the  span- 
wise  macroscale,  A3,  for  both  the  open  and  closed 
channel  cases.  Recall  that  the  solid  wall  is  at  Xt  =  0 
and  the  free  surface  (or  centerline)  is  at  xt  =  1.  The 
scales  are  nominally  given  in  terms  of  outer  units 
(i.e.,  channel  height)  since  conversion  to  wail  vari¬ 
ables  can  easily  be  obtained  by  multipyling  by  13.5  in 
the  open  channel  case  and  125  for  the  closed  chan¬ 
nel.  For  the  streamwise  macroscales,  significant  dif¬ 
ferences  between  these  two  cases  are  apparent  in  the 
rather  large  region  0  4  <  X2  <  1.0.  Here,  the  most 
notable  observation  is  that  the  U3  velocity  compo¬ 
nent  Ai  length  scale  changes  by  a  factor  of  approxi¬ 
mately  three  (as  0  60  at  the  open  channel  free  surface 
relative  to  019  at  the  centerline  of  the  closed  chan¬ 
nel).  For  the  streamwise  velocity  component,  Ai  dif¬ 
fers  only  slightly  at  the  free  surface  from  the  closed 
channel  centerline  value  (1.4  vs.  1.14);  however,  there 
is  a  substantially  different  behavior  in  how  these  final 
values  are  attained.  The  closed  channel  length  scale 
smoothly  achieves  an  asymptotic  value  of  order  one 
at  the  centerline,  whereas  the  behavior  near  the  free 
surface  resembles  that  observed  near  the  solid  bound¬ 
ary  (i.e.,  a  peak  at  some  distance  from  the  boiuid- 
ary  indicative  of  the  source  layer  thickness).  The 
free  surface  effects  on  the  spanwise  macroscales,  .\3, 
are  confined  to  a  smaller  region  (0.8  <  X2  <  10) 
than  those  on  the  streamwise  macroscales.  Substan¬ 
tial  differences  are  again  observed  between  the  values 
attained  at  the  free  surface  relative  to  those  at  the 
closed  channel  centerline.  Both  the  ui  and  U2  compo¬ 
nent  macroscales  differ  by  a  factf-r  of  about  two  with 
the  ui  scale  larger  and  the  »  scale  correspondingly 
smaller  at  the  free  surface. 
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Ai;  =  J  Rjj{xi)dxi.  (8) 

0 

It  should  be  noted  that  in  some  circumstances,  par¬ 
ticularly  for  the  streamwise  velocity  component,  the 
correlation  function  Rjj{Axi)  does  not  decay  suffi¬ 
ciently  at  the  end  of  the  computational  domain  so 
that  the  macroscale  given  by  (8)  may  underestimate 
the  true  eddy  length  scale.  Also,  since  quasi-periodic 
structures  exist  close  to  the  wall  with  their  periodic¬ 
ity  primarily  in  the  spanwise  X3  direction,  RjjiAxj) 
can  be  negative.  This  also  has  the  effect  of  produc¬ 
ing  a  macroscale  which  underestimates  the  length  of 


In  Fig.  5  the  microscale  results  are  compared 
for  the  open  and  closed  channel  flows.  Trends  similar 
to  those  observed  for  the  macroscale  results  of  Fig. 
4  are  apparent.  ”  he  streamwise  scale,  A],  of  the  «3 
component  is  'i:ger  at  the  free  surface  than  at  the 
closed  channel  centerline  (0.40  vs.  0.29):  whereas. 
A]  for  the  u;  component  is  somewhat  lower  (0  21  vs 
0.27).  Similar  to  the  spanwise  macroscale  results,  the 
spanwise  microscales  differ  at  the  free  surface  with 
the  uj  component  scale  larger  and  the  «2  component 
scale  smaller  by  the  same  amounts  than  their  clo.sed 
channel  centerline  values. 

The  vertical  macroscale,  .V;;.  is  given  by; 
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A2;(r2)=  r  R„i0,0,X2.x'^)dx'.,.  (9) 

Jo 

At  a  given  vertical  location  x^,  A2  gives  a  measure 
of  the  vertical  size  of  a  typical  eddy  that  can  ex¬ 
ist  at  that  depth.  These  results  are  shown  in  Fig. 
6.  It  is  observed  that  A21  and  A22  increase  contin¬ 
uously  away  from  the  solid  wall  until  approximately 
X2  =  0.6,  where  the  free  surface  effects  become  evi¬ 
dent  as  both  decrease.  Generali)  speaking,  both  the 
free  surface  and  solid  boundaries  have  similar  overall 
effects  in  that  the  eddy  size  that  can  exist  there  is 
smaller  than  that  existing  away  from  the  boundary. 
However,  the  lack  of  viscous  dissipation  near  the  free 
surface  apparently  allows  for  a  somewhat  greater  ver¬ 
tical  extent  than  that  found  near  a  solid  wall,  as  one 
would  expect.  The  A23  macroscale  results  may  not 
be  easily  interpreted  near  the  solid  boundary  since 
Rsoixi.  ■rz')  becomes  negative  there,  presumably  due 
to  the  presence  of  the  counter-rotating  vortex  struc¬ 
ture  typically  associated  with  the  wall  layer.  Near 
the  free  surface  however,  Rooixi.x^')  is  strictly  posi¬ 
tive  so  that  the  interpretation  of  A23  as  a  length  scale 
is  more  meaningful.  As  with  the  other  two  vertical 
macroscale  components,  the  proximity  of  the  free  sur¬ 
face  (or  perhaps  boundary  in  general)  evidences  itself 
as  a  decrease  in  the  A23  scale.  These  results  generally 
tend  to  confirm  the  HG  model  prediction  of  a  strong 
truncation  of  the  vertical  extent  of  a  typical  eddy  near 
the  free  surface. 

To  summarize,  the  macroscale  results  are  partic¬ 
ularly  useful  in  understanding  the  change  in  shape 
of  a  typical  eddy  as  it  interacts  with  the  free  sur¬ 
face.  Predominantly,  the  free  surface  effects  relative 
to  the  closed  channel  centerline  behavior  are  larger 
spanwise  scales  associated  with  the  streamwise  veloc¬ 
ity  component  ui;  generally  smaller  streamwise  and 
spanwise  scales  for  the  wall  normal  velocity  compo¬ 
nent  ^2;  and  larger  streamwise  scales  associated  with 
the  spanwise  velocity  component  U3.  Additionally, 
the  vertical  macroscales,  A2;,  associated  with  all  three 
velocity  components  decrease  on  approach  to  the  free 
surface  This  presents  a  reasonably  clear  picture  of 
eddies  which  flatten  out  or  become  pancake-like  in 
the  horizontal  plane  as  they  undergo  a  reduction  in 
their  vertical  extent.  This  is  certainly  in  agreement 
with  one's  intuitive  expectation  of  the  similar  effect  of 
a  solid  boundary  on  an  impinging  eddy.  These  results 
indicate  a  streamwise  stretching  by  a  factor  of  about 
three  and  aspanwiseelongationof  approximately  two. 

It  does  not  seem  possible  to  make  a  quantita¬ 
tive  comparison  of  these  results  with  the  HG  model 
in  its  present  form,  since  their  model  does  not  ac¬ 
count  for  the  strong  anisotropy  of  the  open  channel 
flow  studied  here.  Obviously  there  is  some  difficulty 
in  defining  a  far-field  integral  length  scale  in  tliis  flow. 


whereas  this  scale  is  well  defined  for  homogeneous 
turbulence.  However,  since  the  source  layer  of  the 
HG  model  is  generated  by  the  no  mass-flux  boundary 
condition,  it  would  seem  that  the  scales  of  U2  should 
give  the  best  indication  of  the  extent  of  free  surface 
effects.  In  Fig.  7  the  mean-square  turbulence  inten¬ 
sities  scaled  by  the  local  turbulent  kinetic  energy  are 
shown.  This  demonstrates  clearly  that  all  three  veloc¬ 
ity  components  show  the  effects  of  surface  proximity 
at  a  distance  of  about  Ar2  =  0.3.  (See  Swean  et 
for  a  discussion  of  the  redistribtion  of  energy  from 
the  vertical  velocity  component  to  both  the  stream- 
wise  and  spanwise  components.)  Indeed  the  vertical 
velocity  component  is  the  only  one  which  has  asso¬ 
ciated  macroscales  of  about  this  size  and  which  are 
smaller  than  the  macroscales  at  the  centerline  of  the 
closed  channel  flow.  This  seems  to  indicate  that  the 
length  scales  for  U2  also  determine  the  extent  of  the 
source  layer.  This  result  is  certainly  consistent  with 
the  source  layer  model  of  HG  as  noted  above. 

5.  WALL-LAYER  STREAK  SPACING 

In  Section  3  the  persistence  of  a  spanwise  peri¬ 
odic  structure  at  relatively  large  distance  away  from 
the  wall  was  noted  for  the  open  channel  flow.  These 
streamwise-elongated  structures,  commonly  referred 
to  as  wall-layer  streaks,  appear  in  flow  visualization 
studies  as  regions  of  low-speed  fluid  close  to  wall. 
Though  there  remains  some  controversy  as  to  the  sig¬ 
nificance  of  the  streaks,  there  appears  to  be  increasing 
evidence  that  they  are  indicators  of  quasi-streamwise 
oriented  vortices.  These  vortices  in  turn  are  thought 
to  play  a  role  in  the  production  of  new  turbulence 
and  in  Reynolds  stress  production.  They  were  first 
observed  experimentally  by  Hama  (see  Corrsin*^)  and 
later  studied  in  more  detail  by  Kline  el  alJ*.  Their 
visualization  studies  showed  that  the  streaks  were 
typically  observed  below  xj  =  30  and  that  they  oc¬ 
curred  randomly  in  space  and  time.  Kline  ei  al.  also 
found  that  the  average  spanwise  spacing  between  the 
streaks.  A'*’ ,  was  approximately  100  essentially  in¬ 
dependent  of  Reynolds  number.  The  experimental 
results  reported  by  Nakagawa  and  Nezu'^  for  open 
channel  flow  indicate  that  the  mean  streak  spacing 
increases  with  distance  from  the  wall  and  ultimately 
approaches  a  value  of  A+  ss  2X2  for  Xj  >,50.  This  led 
them  to  speculate  that  the  increase  in  length  scale 
resulted  from  a  coalescence  (similar  to  the  pairing  in¬ 
teraction  observed  in  free  shear  flows)  of  neighboring 
low-speed  streaks  as  the  distance  from  the  wall  in¬ 
creases.  However,  it  should  be  noted  that  for  the 
locations  above  X2  >  30  the  spanwise  length  scale 
they  observe  is  very  weak  and  may  not  necessarily 
correspond  to  well-defined  streaks. 

More  recently,  Smith  and  Metzler’®,  in  agree¬ 
ment  with  the  findings  of  Nakagawa  and  Nezu,  found 
that  the  average  spanwise  wavelength  increased  from 
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93  at  X2  =  ^  *2  =  30.  They  noted,  how¬ 

ever,  that  beyond  z}  =  30  that  the  streaks  were 
not  sufficiently  well  defined  to  warrant  making  streak 
counts  and  they  suggest  that  Xj  w  40  is  the  upper 
limit  for  which  extended  regions  of  low-speed  fluid 
continue  to  exist.  Also  in  agreement  with  Nakagawa 
and  Nezu,  they  found  that  streak  coalescence  in  the 
region  10  <  xt  <  30  contributed  to  the  increasingly 
disrupted  streak  pattern  and  overall  increase  in  the 
spanwise  length  scale  with  distance  from  the  wall. 
They  noted  that  since  the  most  active  merging  occurs 
in  the  region  of  maximum  turbulent  energy  produc¬ 
tion,  this  merging  process  may  very  well  be  important 
to  the  turbulence  production  cycle.  The  observation 
of  wall-layer  streaks  has  not  been  confined  strictly  to 
experimental  studies;  various  numerical  simulations 
appear  to  very  satisfactorily  capture  the  wall-layer 
dynamics.  As  an  example,  K.MM  show  very  good 
agreement  with  the  experimental  determinations  of 
the  variation  of  mean  spanwise  streak  spacing  with 
distance  from  the  wall. 

In  Fig.  8  the  dependence  of  the  mean  streak 
spacing  on  xj  is  presented  for  both  open  and  closed 
channel  flows  along  with  various  experimental  and 
numerical  simulation  results.  Here  the  mean  streak 
spacing  is  defined  as  twice  the  spanwise  distance  at 
which  f?ii(Ax3)  reaches  a  minimum  (note  that  this 
definition  of  A+  is  consistent  with  KMM).  The  results 
indicate  clearly  that  for  xj  <  12  there  is  excellent 
agreement  among  all  the  studies  that  A+  ~  100.  The 
streak  spacing  in  the  open  channel  case  shows  a  jump 
from  a  value  of  105  at  xj  =  12  to  approximately  130 
at  Xj  =  15.  Farther  from  the  wall,  the  open  chan¬ 
nel  streak  spacing  increeises  at  a  rate  which  is  roughly 
the  same  as  in  the  closed  channel  but  always  remains 
larger.  Thus,  Fig.  8  shows  that  at  x^  —  65,  A+  has 
attained  a  value  roughly  twice  twice  that  observed  at 
x^  ~  12.  It  is  also  clear  that,  unlike  the  closed  chan¬ 
nel  case  where  the  streak  spacing  is  apparently  only 
clearly  defined  out  to  perhaps  xj  ~  50,  the  spanwise 
periodic  streaks  can  definitely  be  observed  well  out 
to  x^  ~  80  in  the  open  channel  flow.  This  is  in  fact 
within  the  logarithmic  layer  of  the  streamwi.se  veloc¬ 
ity  profile  (see  Swean  et  al}°  for  the  open  channel 
flow  mean  velocity  profile).  It  is  apparent  that  the 
the  streaky  structure  in  open  channel  turbulence  is 
both  larger  in  scale  and  persists  farther  from  the  wall 
than  in  closed  channel  turbulence. 

Though  the  reasons  for  these  differences  are  far 
from  evident,  a  few  speculations  are  now  offered  to 
explain  this  behavior.  The  only  obvion.s  difference  be¬ 
tween  these  two  flow  fields  is  the  boundary  condition 
imposed  on  the  upper  boundary.  Since  tlie  shear-free 
boundary  suppresses  the  production  of  new  turbu¬ 
lence  we  can  envision  that  this  must  in  turn  suppress 
the  ejection  of  low  momentum  fluid  that  would  oth¬ 


erwise  occur  if  the  boundary  were  rigid  (i.e  no-slip 
boundary  conditions).  The  suppression  of  these  ejec¬ 
tions  from  the  top  boundary  (free  surface)  toward  the 
bottom  boundary  (solid  wall)  may  be  responsible  for 
slowing  down  the  production  of  turbulence  by  low¬ 
ering  the  probability  of  shear  layer  formation  in  the 
region  of  the  solid  wall.  Thus,  the  larger,  more  coher¬ 
ent  streaks  near  the  wall  in  open  channel  turbulence 
ate  due  to  the  suppression  of  turbulence  production 
at  the  free  surface.  Furthermore,  in  boundary  layer 
flows  the  entrainment  of  outer  irrotaional  fluid  may 
act  like  ejections  from  the  upper  wall  of  a  channel 
In  this  sense,  boundary  layer  flow  has  a  closer  resem¬ 
blance  to  closed  channel  flow  than  to  open  channel 
flow.  These  results  suggest  that  the  outer  flow  does 
have  an  effect  on  the  wall  region  at  least  at  these  low 
Reynolds  numbers. 

6.  CONCLUSIONS 

The  structure  of  turbulence  near  a  free  surface 
has  been  studied  using  the  results  of  a  direct  simula¬ 
tion.  The  method  employed  here  to  understand  this 
structure  is  to  compare  the  open  channel  turbulence 
results  with  its  well  studied  closed  channel  turbu¬ 
lence  counterpart.  A  comparative  examination  of  the 
turbulent  macroscales  and  microscales  in  these  two 
cases  reveals  a  significant  flattening  of  a  typical  eddy 
near  the  free  surface.  This  flattening  is  evidenced 
by  a  notably  larger  streamwise  scale  associated  with 
the  spanwise  velocity  component  and  a  comensurately 
larger  spanwise  scale  of  the  streamwise  velocity  com¬ 
ponent.  Additionally,  the  vertical  macroscales  for  all 
three  velocity  components  are  smaller  than  the  com¬ 
panion  closed  channel  values.  In  this  flow,  the  source 
layer  described  by  the  Hunt-Graham  model  appears 
to  extend  about  0.3  channel  heights  below  the  free 
surface.  This  source  region  seems  to  correlate  most 
strongly  with  the  structure  of  the  vertical  component 
of  velocity.  All  length  scales  for  this  component  are 
about  the  correct  size  and  all  are  smaller  than  the 
macroscales  at  the  centerline  of  the  closed  channel 
flow.  A  direct  quantitative  comparison  of  these  re¬ 
sults  with  the  HG  model  is  not  possible  since  the 
model  does  account  for  the  strong  anisotropy  of  the 
turbulence  present  in  this  flow.  The  energy  spectra 
for  all  three  components  of  velocity  are  in  qualita¬ 
tive  agreement  with  the  HG  model  though  resolution 
effects  and  anisotropy  limit  direct  quantitative  com¬ 
parison. 

An  interesting  phenomenon  revealed  by  the  cur¬ 
rent  study  is  the  increase  in  size  and  persistence  of 
the  spanwise  periodic  structure  near  the  wall.  In  open 
channel  turbulence  this  structure  is  larger  in  scale  and 
penetrates  farther  into  the  flow  than  in  closed  channel 
turbulent  flow.  The  origin  of  this  effect  is  not  clear 
but  one  possible  explanation  is  that  the  free  surface 
suppresses  interactions  which  would  normally  occur 
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between  the  no-slip  boundaries  of  channel  flow.  This 
observation  lends  support  to  the  possibility  that  the 
outer  flow  has  a  significant  effect  on  wall  layer  struc¬ 
ture.  In  future  work,  a  quantitative  comparison  of 
these  results  with  a  modified  form  of  the  HG  model 
will  be  attempted  and  higher  resolution  simulations 
will  be  undertaken  to  further  elucidate  the  structure 
of  the  turbulence  near  the  free  surface. 
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Fig,  1  Streamwise  two-poiiit 
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Fig.  8  Variation  of  mean  spanwise  streak  spacing  with  wall  normal  distance.  Present  study:  x ,  open  channel, 
R'  =  135;  +,  closed  channel,  R'  =  125.  Smith  and  Metzler*®:  o,  boundary  layer  expt.,  Rg  =  2030.  Schraub 
and  Kline'^:  □,  boundary  layer  e.xpt.,  Re  =  1325.  Nakagawaand  Nezu'^:  open  channel  expt.,  y,  R'  =  248; 
A,  R‘  =  378.  Kim,  Moin,  and  oser'’:  o,  closed  channel  simulation,  R'  =  180. 


